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Abstract. We have analyzed the galaxy number density 
and luminosity density profiles of rich clusters of galaxies 
from redshifts z^O to z~0.5. We show that the luminosity 
profile computed with bright galaxies (Mr < —21) is sig- 
nificantly cusped in the center of the clusters, whatever the 
redshift. This is in agreement with the dark matter pro- 
files predicted by numerical simulations. The galaxy num- 
ber density profile for the bright galaxies is fitted equally 
well with a core model or a cusped model. In contrast, 
the luminosity and the galaxy number density profiles of 
the fainter galaxies are significantly better fitted by a core 
rather than a cusp model. We did not detect any statisti- 
cally significant different fits when applied to data in range 
from z^O to z~0.5. The difference in profile between faint 
and bright galaxies may be due to the rapid (relative to 
the age of the universe at z—0 versus z—0.5) destruction of 
the faint galaxies by tidal forces and merging events in the 
denser central regions of the clusters. This process could 
erase the cusp by turning faint galaxies into diffuse light. 
In this case, the galaxies (with a cusp visible in the bright 
galaxy number density and mainly in luminosity profiles) 
would trace the total mass distribution. 
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1. Introduction 

Clusters of galaxies are key cosmological probes. They 
are one of the main products of the hierarchical models 
and the largest virialized structures in the Universe. Re- 
cent studies (e.g. Navarro et al. 1997) argue that whatever 
the cosmological model, a universal profile (NFW profile) 
with a cusp at the center should describe the cluster dark 
matter profile, in contrast to the beta-model (King 1962) 
which has been used from many years (cf Sarazin 1986). 
The shape of the galaxy number and luminosity density 
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profiles (and their dependence on redshift) relates directly 
to the physical processes acting at the center of the clus- 
ters and we propose to investigate in this way these pro- 
cesses. While the dark matter profile (close to the total 
mass profile) seems to be very concentrated (Navarro et 
al. 1997), the galaxy distribution is nearly flat in the cen- 
ter of the nearby clusters (Adami et al. 1998: ENACSVII). 
Furthermore, the X-ray gas distribution cannot be used to 
distinguish between a flat and a cusped model except for 
the nearby clusters (e.g. Durret et al. 1994) due to a lack 
of resolution with Rosat, Einstein and Asca and due to 
cooling-flows. 

In this paper we contribute to the effort to solve this 
puzzle by studying the galaxy number density profile and 
the galaxy luminosity density profile (the sum of the lu- 
minosity inside the galaxies; we will consider the diffuse 
light only at the end of the paper). This approach is clearly 
not new (e.g. Mazure et al. 1986), but we used here very 
large, homegeneously analyzed, samples. The use of all 
the optical wave band galaxy light allows us to take into 
account merging events. This method enables us to re- 
cover " erased" cusps as it is the light of merged galaxies 
which is counted, not the number of galaxies. The best 
way to investigate the mass profile would be to compute 
this profile via the velocity dispersion profile (e.g. Carlberg 
et al. 1997a and b, Biviano et al. 2001), the X-ray tem- 
perature profile, or gravitational lensing, but these tech- 
niques require large amounts of data, and therefore, it is 
very difficult to obtain large samples of clusters. Moreover, 
galaxy luminosity density profiles for early type galaxies 
are claimed (e.g. Kaiser 1999) to be similar to the mass 
profile in clusters. This is an important question as galax- 
ies would be sufficient to trace the mass profile. We will use 
in this paper homogeneous samples of thousands of galax- 
ies and tens of clusters to compare the projected galaxy 
number density and luminosity profiles, and the simulated 
dark matter profiles (e.g. Navarro et al. 1997). 

In Section 2, we describe the samples. In Sections 3 
and 4, we describe the methods we used and the profiles 
we generated. Section 5 discusses the results and Section 
6 gives a summary of the results. 
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To relate values of z to distance and to be in agreement 
with Adami ct al. (1998, ENACSVII), we have taken H 
= 100 km.s^Mpc" 1 , q = and A = 0. 

2. Samples 

2.1. Low redshift sample 

We used the COSMOS composite cluster which is de- 
cribed in detail in ENACSVII, but we briefly enumerate 
the main points here: we compiled photometric COSMOS 
data (Heydon-Dumbleton ct al. 1989) for 77 of the 107 
ENACS clusters (z<0.1). We have adopted a limiting mag- 
nitude of bj — 20. 

We have selected the clusters with: 

1. a low level of substructure (both visually and with a 
Dressler-Shectman test) in a 10 core radii on a side 
square (^1200 kpc h _1 ), with more than 10 known 
rcdshifts in the main group (to avoid superposition ef- 
fects) ; 

2. a redshift lower than 0.1 

3. and a converging solution of the fit of analytical density 
profiles. 

These selection criteria lead to a sample of 29 homo- 
geneously selected clusters. With this sample we built a 
composite cluster with about 5000 galaxies. For each of 
the 29 clusters, we determined the center by a Maximum 
Likelihood fit and we confirmed the results with other es- 
timators (X-Ray centers, cD galaxies, see also Ulmer et 
al. 1992). For each galaxy, the distance from the center 
has been scaled both with the core radius of the cluster 
and the r2oo (radius where the density is equal to 200 
times the critical density) in order to take into account 
that different clusters have different sizes (see ENACS VII 
for a discussion of this point). We have taken the elon- 
gation of the clusters into account by 'circularizing' the 
individual galaxy distributions by increasing all projected 
distances orthogonal to the major axis, thus 'expanding' 
the distribution parallel to the minor axis by a factor de- 
duced from the ellipticities of individual clusters. This cor- 
rection is important because the superposition of galaxy 
distributions with randomly distributed orientations will 
cause outer densities to be underestimated with respect 
to the inner ones, which is producing an artificial cusp (cf 
ENACS VII). 

2.2. High redshift sample 

We used a sample of 7 clusters from the COP (Adami et 
al. 2000a, Holden et al. 2000) and CNOC (e.g. Carlberg 
et al. 1996) surveys selected with the same conditions as 
the low redshift sample: 

1. low level of substructures in a 10 core radii square 
(~1200 kpc) from visual analysis 



2. more than 10 redshifts in the main group 

3. converging solutions of the fit of analytical galaxy 
number and luminosity density profiles (this excludes 
2 clusters: MS1008 and MS1224). 

This sample is described in Table 1: the redshift range of 
this sample of distant clusters is [0.33;0.55] (mean redshift 
= 0.42). We have built a composite cluster of 459 galaxies 
with these 7 clusters, correcting for ellipticity and orien- 
tation of each real cluster, in a similar manner as we did 
for the nearby sample. 

We note that MS 1358 is regarded in the CNOC papers 
as a complex cluster. However, we used in this paper only 
the central regions of this cluster in order to focus on the 
central cluster shape: we limited this cluster to a region 
of ^160" radius. With this restriction, MS1358 can be 
included in our sample. We stress, however, that using 
larger areas will induce significant levels of substructures. 

3. Methods and Results 

3.1. Low redshift galaxy number density and luminosity 
profiles 

The individual parameters (with individual clusters) of 
the low redshift sample profiles are described in Adami 
ct al. (1998: ENACSIV). They have been computed with 
a Maximum Likelihood fit. We can easily generalized the 
same code to fit a luminosity profile for the composite 
clusters. We have weighted each galaxy with its lumi- 
nosity (assuming all the galaxies at the cluster redshift): 
Lunik- Similarly to ENACSVII, the probability that the 
assumed profile 'produces' a galaxy in position (x k ,y k ) is 
o~{x kl y k ) (with k=l...N). The combined probability that 
the assumed profile will produce galaxies in the positions 
( x fc,yfe) that they actually have is: 

N 

L= f] Lum k x a(x k ,y k ) 
fe=i 

The best fit model is found from a maximization of L. 
We have chosen, instead, to compute -ln(L) with: 

N 

In L = J2 ln(Lum k x (j(x k ,y k )) 

k=l 

The best model (which produces the best fit) is found 
from a minimization of -ln(L). We minimize this value by 
using the MINUIT package (e.g. ENACSVII for details). 

To be able to quantify the degree of cuspiness at the 
center of the clusters, wc fitted in ENACSVII both a pro- 
jected beta- model: 

Kr)=ti (l/(l + (?-J 2 r+ti b 

and a cusped profile: 

Mr)= M o(l/(£(l + £) 2 ))' 3 +M 6 

For the luminosity profiles, fi is the surface brightness, 
which replaces the surface density a of the formulae in 
ENACSVII. These models, fitted to the individual clus- 
ters, have 7 free parameters: two parameters for the posi- 
tion of the center (x ,yo), two parameters to describe de- 
viations from symmetry (ellipticity e and position angle <j>, 
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included in r), two parameters that specify the profile (r c 
and (3) and the background density /i& (assumed constant 
within the aperture of each cluster). 

To be able to compare the fit quality of the two models 
(cusped profile and beta- model), we use the logarithm of 
the ratio of the Maximum Likelihood values for each fit as 
explained in ENACSVII. This value has a x 2 distribution 
(Meyer 1975) which we then use to determine the statis- 
ticial significance level of the best fit of one profile type 
versus the other. 

Moreover, we have checked the consistency between 
the center of the clusters determined by using the galaxies 
(see ENACSVII) and the center determined by using the 
galaxies weightened with their luminosity. We find a shift 
of only (32 ± 21) kpc. This shift is not large enough to 
erase a cusp (e.g. ENACSVII Fig. 8) and we have assumed 
the centers using only the galaxy number density. 

3.2. High redshift galaxy number density and luminosity 
profiles 

We used a similar approach as for the low redshift sample, 
but with an additional step to take into account the high 
background level. This is because it is nearly impossible to 
measure redshifts for all galaxies of all the clusters in our 
high redshift sample due to the large amount of observ- 
ing time needed. Therefore we used the red sequences in 
the Color Magnitude Relations (CMR hereafter) to reject 
many of the field galaxies along the lines of sight (e.g. Yee 
et al. 1999). 

The color intervals we used are given in Table 1. 
We chose these intervals manually in order to match the 
prominent structures in the CMR of each line of sight. 
We have confidence that we made good choice of inter- 
vals as these intervals are very similar to the theoretical 
predictions of Kodama & Arimoto (1997). 

We give in Table 1 the best fit values of r c and (3 for 
the 7 high redshift clusters, the 2 analytical models (beta- 
model and cusped profile) and the galaxy number and 
luminosity density profiles parameters (r c is ranging from 
63 to 299 kpc and (3 from 0.8 to 1.04 for a beta-model). 
We give also the fit quality: which analytical model fits 
best, core or cusp? 

We now examine whether that these values could have 
been affected by the following bias: the CMR rejects not 
only field galaxies but also the late type galaxies that are 
cluster members. Since these galaxies are preferentially lo- 
cated in the outskirts of clusters (e.g. Adami et al. 1998a), 
this could affect the density profiles. 

- At low redshift, this is not a serious concern as more 
than 75% of the cluster galaxies in a 1200 /i _1 kpc square 
are early type objects (Adami et al. 1998a). 

- For the higher redshift clusters, however, the spiral 
fraction increases (e.g. Dressier et al. 1997) and can con- 
tribute for up to 60% in the outer parts of the clusters 
(Ellingson et al. 2000). As we want to focus on the central 
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Fig. 1. Significance level of the core/cusp discrimination 
in our simulations (y-axis: percentage) as a function of the 
field contribution to the total number of galaxies along the 
line of sight (x-axis: percentage). The horizontal line is the 
70% significance level. 

cluster shape, we limited ourselves to a region of ~0.5r2oo 
(5 core radii) . In this region, the field- like galaxy contribu- 
tion is lower: ~30%. We have tested, however, the signifi- 
cance of the effect by removing these galaxies when we fit 
the density profiles. We used the CNOC cluster MS1358 
because redshifts are available for nearly all the galaxies 
along the line of sight. We can, therefore, select the cluster 
member galaxies on a firm basis. We see in Table 1 that 
using the CMR or the redshift selection does not affect in 
a significant manner the fitted parameters: these param- 
eters are similar at the 1-cr level: for a beta-model r c is 
ranging from 117 to 95 kpc and (3 from 0.91 to 0.98. This 
is not a definitive test as we are using a single cluster, but 
the results are still very suggestive. 

Using the CMR seems, therefore, a good way to remove 
field galaxies in order to discriminate between cusped 
profiles and profiles with a core. However, as we go to 
faint magnitudes, the uniform field contribution becomes 
stronger and stronger. This field contribution is between 
50% and 75% for the clusters in our sample (estimated 
from the redshift catalogs of these cluster lines of sight). 
This could have the effect that the cluster density profiles 
will appear flatter for faint galaxies. 

-This does not affect our results for nearby clusters 
(we do not sample the clusters at magnitudes fainter than 
R=19, where the CMR becomes to be less contrasted com- 
pared to the field contribution: Adami et al. 2000b). 

-This does not affect our results for distant clusters 
(z^0.4) for the brightest bin of Table 3 (magnitudes 
brighter than R=19). 

-This could have an influence for the 3 faintest bins 
of the distant composite cluster. However, the goal of this 
paper is to discriminate between density profiles with a 
cusp or with a core. In order to quantify the probabil- 
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Table 1. High rcdshift sample with cluster name, selection of the cluster galaxies (color interval or redshift), redshift, 
characteristic radius and (3 for the beta-model (galaxy number density profile: dens, gal.), characteristic radius and 
(3 for the cusped model (galaxy number density profile: dens, gal.), characteristic radius and (3 for the beta-model 
(galaxy luminosity profile: dens, lum.), characteristic radius and (3 for the cusped model (galaxy luminosity profile: 
dens, lum.) and best fit (galaxy number and luminosity density profiles). The quoted errors are at the 1-cr level. No 
error means that we were not able to get a reliable estimate for this error. 



Name 


Selection 


z 






r c dens. gal. beta-model 


(3 dens. gal. beta-model 


r c dens. gal. cusped 


j3 dens. gal. cusped 


best fit dens. gal. 


r c dens. lum. beta-model 


(3 dens. lum. beta-model 


r c dens. lum. cusped 


(3 dens. lum. cusped 


best fit dens. lum. 


MS0302 


g-R: [1.00,1.80] 


0.42 






63±12 kpc 


0.93±0.06 


93±32 kpc 


0.52+0.04 


beta-model: not significant 


62±12 kpc 


0.94±0.09 


94±29 kpc 


0.52+0.05 


beta-model: not significant 


MS0451 


g-R: [1.25,1.85] 


0.54 






121±32 kpc 


0.87±0.17 


209±108 kpc 


0.53+0.07 


beta-model: not significant 


1 1 7+29 kne 


88+0 16 


210+120 kne 


53+0 08 

kj • i yj • \jyj> 


hptfl-tnnrlpl • not ^ioTiifipflnt 


MS1358 


redshift 


0.33 






117±26 kpc 


0.98±0.09 


455±97 kpc 


0.81+0.05 


cusped: not significant 


119±29 kpc 


0.99±0.08 


469±95 kpc 


0.81+0.05 


cusped: not significant 


MS1358 


g-R: [1.02,1.50] 


0.33 






95±22 kpc 


0.91±0.07 


516+172 kpc 


0.83+0.06 


cusped: not significant 


101±13 kpc 


0.92±0.04 


509+151 kpc 


0.83+0.05 


cusped: not significant 


MS1621 


g-R: [0.85,1.85] 


0.43 






67 kpc 


0.80 


735 kpc 


0.86 


beta-model: not significant 


68 kpc 


0.81 


737 kpc 


0.85 


beta-model: not significant 


PDCS16 


V-I: [1.35,2.11] 


0.40 






299±39 kpc 


0.99±0.20 


882 kpc 


0.68 


beta-model: not significant 


303±36 kpc 


0.98±0.18 


885 kpc 


0.67 


beta-model: not significant 


PDCS30 


V-I: [1.10,1.90] 


0.33 






189 kpc 


1.04 


302+67 kpc 


0.53+0.04 


beta-model: 85% level 


189 kpc 


1.01 


227+69 kpc 


0.46+0.03 


beta-model: 85% level 


PDCS62 


V-I: [1.24,2.20] 


0.46 






196±10 kpc 


1.01±0.26 


284+73 kpc 


0.53+0.04 


beta-model: 85% level 


192 kpc 


1.02 


282+71 kpc 


0.53+0.03 


beta-model: 85% level 



ity to have a cusped profile mis-interpreted as a profile 
with a core due to a high field contribution, we performed 
simulations. We generated a circular artificial cluster pro- 
file with a cusp in the center (as described in Section 3.1 
with r c =200kpc and (3=0.67) and we added to this cluster 
a uniform field contribution with various galaxy densities. 
We computed finally the relative fit quality between a pro- 
file with a core and with a cusp. For field contributions as 
high as 80% of the total number of galaxies along the line 
of sight, a model with a cusp is prefered at the 80% sig- 
nificance level. For contributions equal to 85%, a model 
with a cusp is prefered at the 70% significance level. For 
higher field contributions, we are not able to discriminate 
between a profile with a core or with a cusp (a model with 
a cusp is still prefered but with significance levels lower 
than 65%, which is not significant). The results are plotted 
in Fig. 1. This means that even with a high field contri- 
bution, we are still able to discriminate between profiles 
with a core or with a cusp because the field contribution 



along the lines of sight we used is not higher than 80%. 
Our results are, therefore, valid regarding the core/cusp 
discrimination. 

3.3. Analysis of the results 

-We compare the characteristic parameters r c and [3 of 
the distant clusters with those of the nearby sample 
(ENACSVII). The results are summarized in Table 2 and 
we give the description of these results: we used first the 
individual clusters. We do not see a statistically significant 
evolution with redshift: the agreement between z~0 and 
z^0.4 fits is consistent at the 2a level (for a beta-model, in- 
dividual clusters and galaxy number density, r c ~ ISOkpc 
and (3 ~ 1). If we compare now the fitted parameters r c 
and (3 of the galaxy number and luminosity density pro- 
files (regardless of the redshift), they also agree at the 2-a 
level (for a beta-model, composite clusters and luminosity 
profiles, r c ~ llOkpc and (3 ~ 1.05). 
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Table 2. Galaxy number density and luminosity profile parameters (characteristic radius r c and slope 0) for nearby 
and distant clusters. The values for the individual cluster galaxy number density profiles are the mean and the 1-er 
uncertainty between the individual fits. The values for the composite clusters are the fitted values. No error means 
that we were not able to get a reliable estimate for this error. 



DENSITY 


r c beta-model 


r c cusped 


/3 beta-model 


(3 cusped 


best fit 


Nearby sample (z~0.07) 
(Individual clusters) 


<128>±88 kpc 


<292>±191 kpc 


<1.02>±0.08 


<0.61>±0.05 


beta-model 
Not significant 


Distant sample (z~0.4) 
(Individual clusters) 


<147>±79 kpc 


<431>±268 kpc 


<0.94>±0.08 


<0.64>±0.14 


beta-model 
Not significant 


Nearby sample (z~0.07) 
(Composite cluster) 


89±5 kpc 


318±34 kpc 


1.00±0.02 


0.56±0.01 


beta-model 
99% level 


Distant sample (z~0.4) 
(Composite cluster) 


119±8 kpc 


215±9 kpc 


1.08±0.05 


0.65±0.04 


beta-model 
85% level 


LUMINOSITY 


r c beta-model 


r c cusped 


(3 beta-model 


P cusped 


best fit 


Nearby sample (z~0.07) 


101±7 kpc 


281±41 kpc 


1.01±0.03 


0.58±0.03 


beta-model: not significant 


Distant sample (z~0.4) 


119±9 kpc 


232 kpc 


1.08±0.04 


0.66 


beta-model: 85% level 



-As for the nearby sample (e.g. ENACSVII), the clus- 
ter galaxy distribution of the distant cluster sample is bet- 
ter fitted with beta- models than with cusped profiles. Al- 
though the fit is better, it is not significantly better than 
a 75% significance level. 

-To improve the statistical significance of the fits in or- 
der to distinguish between cusped profile and beta-profile 
fits, we used also composite clusters. The method used to 
build these composite clusters is described in Section 2.1. 
We found that the behavior of the fits was similar for the 
distant and nearby samples: the bright galaxies better fol- 
lowed a cusped luminosity profile, while the faint galaxy 
luminosity and galaxy number density profiles were bet- 
ter fitted with beta-models. This is also illustrated for the 
nearby sample with Fig. 2. We see that the bright galaxy 
luminosity density profile seem to be more peaked than 
for the fainter galaxies. 

4. Cores or Cusps: Comparison with Literature 

We have analyzed in this paper a very large sample of 
clusters and we have shown that we have a better fit for 
the galaxy number density if we used a model with a core 
rather than with a cusp (except for the bright galaxies). 
This is apparently contrary to the study of Carlbcrg ct al. 
(1997a) which favored a model with a cusp for the galax- 
ies of the CNOC clusters. We have shown in ENACSVII, 
however, that the results of Carlberg et al. are probably 
explained because Carlberg et al. did not correct for ellip- 
ticity of the cluster profiles, and more importantly, they 
only considered bright galaxies. 

We have shown that if we limit ourselves to bright 
galaxies only (see Tab. 3), there is no statistically signif- 
icant difference between a model with a cusp or with a 
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Fig. 2. Illustration of the luminosity profiles for the r c 
scaled nearby composite clusters (assuming r c ~100 kpc) 
and the four luminosity ranges. The r c axis is in logarith- 
mic units of r c . Each bin has the same number of galaxies 
and the y values are rescaled to the same integral. The 
3 symbols represent the value of the central bin with the 
diffuse light correction applied to the luminosity density 
profiles listed in the figure next to the relevant symbol. 



core for the galaxy number density. This result is in good 
agreement, for example, with the work of Biviano et al. 
(1996) on the Coma cluster. 

For another comparison we remark that Durret et al. 
(1994) showed that, using the X-Ray surface brightness 
profile of 12 clusters (based on ROSAT data), a model 
with a core does not fit significantly better the observa- 
tions than a model with a cusp. 
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Table 3. Best fit for the different models and composite clusters. 



Nearby composite cluster 



Luminosity density profile 



Galaxy number density profile 



All galaxies, r c scaled 
All galaxies, r2oo scaled 
bj [-22.4;-18.78], r c scaled 
bj [-18.77;-18.03], r c scaled 
bj [-18.02;-17.48], r c scaled 
bj [-17.47;-16.89], r c scaled 



beta-model better than cusp: not significant 
beta-model better than cusp: not significant 
cusp better than beta-model 95% level 
beta-model better than cusp 95% level 
beta-model better than cusp 75% level 
beta-model better than cusp 75% level 



beta-model better than cusp 99% level 
beta-model better than cusp 99% level 
beta-model better than cusp: not significant 
beta-model better than cusp 95% level 
beta-model better than cusp 95% level 
beta-model better than cusp 75% level 



Distant composite cluster 



Luminosity density profile 



Galaxy number density profile 



All 


galaxies, 


r c 


scaled 


beta-model better than cusp 85% level 


beta-model better than cusp 85% level 


R 


-22.;-2L], 


r c 


scaled 


cusp better than beta-model 95% level 


beta-model better than cusp: not significant 


R[- 


21. ;-20.5], 


r c 


scaled 


beta-model better than cusp 95% level 


beta-model better than cusp 95% level 


R[- 


20.5;-20.j, 


r c 


scaled 


cusp better than beta-model: not significant 


cusp better than beta-model: not significant 


R 


-20.;-19.], 


r c 


scaled 


beta-model better than cusp 70% level 


beta-model better than cusp 70% level 



7. a - 




7. 2 - 

[-18. 02: -17. 18] 

[-17.47:-16.891 

7 I ' ' ' 1 ' ' ' ' 1 ' ' 1 

-0.5 

Log c rad i lis ) 

Fig. 3. Log. of the ratio of the luminosity profiles and the 
galaxy number density profiles for the r c scaled nearby 
composite clusters (assuming r c ~100 kpc) and the four 
magnitude ranges. 

If we fit the luminosity density profiles for all galaxies 
with a beta-model and a cusped model, the core model is 
still preferred, but the difference is not significant over the 
rcdshift range from to 0.5. For a comparison, we have 
plotted in Fig. 3 the ratio of the galaxy luminosity density 
to the galaxy number density (for the 4 luminosity bins). 
We see that the ratio is constant for the faintest bins, while 
we have an increasing ratio for the brightest bin close to 
the center, inside ~50 kpc. As we have noted in previous 
sections, there seems to be a difference between the distri- 
bution of bright and faint galaxies and also between galaxy 
number density profile and the luminosity density profile. 
It seems that two different shapes apply. The faint galax- 
ies exhibit a core. They have a distribution in agreement 
with a beta-model. The bright galaxies (at least for the 
luminosity profile) show a cusp. This would be in agree- 



ment with the dark matter profiles of the simulations of 
Navarro ct al. (1997). As suggested by Kaiser (1999), the 
bright galaxies would, therefore, trace the cluster mass. 

5. Discussion 

We have presented different arguments favoring either a 
core or cusp in the last sections. The galaxy number den- 
sity profiles do not exhibit a cusp, and the bright galaxies 
cannot be used to distinguish between core and cusp mod- 
els. On the contrary, the luminosity profiles can be used 
to produce a cusp or a core depending on whether just 
bright or faint galaxies are used to generate the luminos- 
ity density profile. This leads to suggest different scenarios 
for the evolution of the bright and faint objects: 

- the bright galaxies luminosity profile could be cusped 
at the formation epoch of the clusters. It could also be- 
come cusped via evolutionary processes, such as a segrega- 
tion process (bright galaxies in the center of the clusters). 

- the faint galaxies, which exhibit a core, could origi- 
nally have been in a cusped distribution, but this distribu- 
tion erased by environmental effects like tidal disruption 
or merging events. Galaxies near the cluster center and on 
radial orbits with low angular momentum would have to 
pass right through the giant ellipticals and would be even- 
tually swallowed up or disrupted in the process. This can 
lead to a lower density of faint galaxies near the cluster 
center. 

5.1. Segregation effects. 

We examine if the luminosity profile cusp for bright galax- 
ies is due to a segregation effect. In Adami et al. (1998a) 
it was shown that the elliptical galaxies are the brightest 
objects in a cluster. Moreover, they are more concentrated 
in the cluster centers than the other morphological types. 
This explains probably why we are not able to distinguish 
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between a core and a cusp when we examine the galaxy 
number density profile of the bright galaxies (ENACSVII) : 
bright galaxies are mainly elliptical and elliptical galaxies 
are mainly in the cluster centers. 

To confirm this effect, we have removed the galax- 
ies brighter than bj = —21 from the brightest bin of 
the nearby sample (28 galaxies). This removes the cD- 
likc galaxies. When re-doing the fit of a beta-model and a 
cusped profile for the bright galaxy bin with the cDs re- 
moved, we still prefer a model with a cusp, with, however, 
at a lower significance level of 90 % (instead of 95 %). This 
effect does not appear, however, to be dominant. 

We now discuss the difference between the galaxy num- 
ber density and the luminosity profiles. A cusp in the 
galaxy number density profile could be erased by merging 
events, leading to a core. However, the luminosity profile 
of the bright galaxies has a cusp, but the same profile for 
the faint galaxies has a core. We discuss a scenario that 
would explain this result in the next subsection. 

5. 2. Destruction of the faint galaxies in the center of the 
clusters? 

Suppose that the formation mechanism of the clusters 
does not initially produce a core, but only a cusp for 
both bright and faint galaxies, as in the simulations (e.g. 
Navarro et al. 1997). This cusp could be "erased" in the 
galaxy number density profile by merging events and " re- 
vealed" in the luminosity profile for the bright galaxies. 
We must now explain why the luminosity profile for the 
faint galaxies has no cusp. This could be explained if we 
consider the other sources of galaxy destruction, beside 
the merging events: the tidal disruptions. Tidal forces in 
the cluster center disrupt part of the faint galaxies, turn- 
ing them into diffuse light, while the bright ones (more 
massive and more robust) are conserved. This explanation 
was, for example, proposed by Merrit (1984). By using 
simulations of the dynamical evolution of the cluster core, 
he predicted a tidal radius (lowest value for the size of a 
not tidally disrupted galaxy) of about 15 h^ 1 kpc. The 
bright galaxies, like cD galaxies, are significantly larger. 
The cusp could be erased in this way only for the faint 
galaxies which are generally smaller. More recent simula- 
tions by Moore et al. (1998) show the same trend for the 
dwarf galaxies to be disrupted in the center of the clus- 
ters. Such a destruction of the faint galaxies has been also 
proposed on observational bases for example by Seeker et 
al. (1997) or Gregg & West (1998). These last authors, 
using deep photometry of the Coma cluster, propose the 
disruption of the faint galaxies as an explanation of the 
lack of dwarf galaxies in the core of this cluster (e.g. also 
Adami et al. 2000b). 

A way to investigate this possibility is, for example, to 
search for cusped profile in the diffuse light. Gregg & West 
(1998) give an approximation of the total luminosity lost 
for the galaxies by disruption and turned into diffuse light 



in the center of the Coma cluster. This value is about 20 
% of the luminosity of a cD galaxy. This contribution is 
around 50 % of the luminosity sum of all galaxies fainter 
than —18.77 in the central bin (where the cusp is sig- 
nificant for the bright galaxies) of our nearby composite 
cluster. After rescaling to the total luminosity to get units 
coherent with Fig. 2, we see on this figure that such a 
contribution is large enough to enhance the three faint 
luminosity profiles, consistently with a cusped model. 

A more quantitative study has been done for this same 
cluster by Bernstein et al. (1995). Using very deep images 
of a cluster-central area of 7'x7', they conclude that the 
faint galaxy luminosity profile is flat or decreasing in the 
central 40 kpc (see also Adami et al. 2000b for a spectro- 
scopic survey of this area), while the luminosity profile of 
the brighter galaxies is peaked. Moreover, considering now 
the diffuse luminosity profile (not associated with visible 
galaxies), they also found a peaked profile. If we argue that 
the cusp of the faint galaxy luminosity profile is erased by 
tidal disruption of part of these galaxies, we should recover 
at least partially this shape in the diffuse light profile. It 
is exactly what Bernstein et al. saw in the Coma cluster. 

5.3. Destruction time scale? 

We demonstrated that the cluster profiles have the same 
central shape from z~0 to z~0.5. This would imply that 
the time scale of the destruction processes described above 
(merging + disruption) would be relatively short. This is 
because all the faint galaxies accreted between z~0 and 
z~0.5 do not contribute to the reconstruction of a cusp. 
The in-falling galaxies seem to be continuously destroyed, 
in a time close to the virialization time which is approxi- 
mately 1.5 Gyear (e.g. Sarazin, 1986) for faint galaxies. 
This estimate is, for example, in good agreement with 
the estimates of the time needed to disrupt a low surface 
brightness galaxy in the simulations of Calcaneo-Roldan et 
al. (2000). An alternative explanation in agreement with 
the observations is that the accretion rate could be already 
very low at z=0.5 (e.g. Ellingson et al. 2000). 

6. Summary 

We have shown in this article that, on one hand, the galaxy 
luminosity density profile is cusped for the bright galaxies 
from z^0 to z~0.5 and exhibits a core for the faint galax- 
ies. On another hand, the galaxy number density profile 
has a core (not significant for the bright galaxies). This 
could be understood if we assume a cusped profile for all 
the galaxies in agreement for example with the simula- 
tions of Navarro et al. (1997). The cusp could be erased 
for the galaxy number density profiles by merging effects 
and for the luminosity profile of the faint galaxies with the 
destruction of these galaxies by tidal forces. This seems to 
be confirmed by two observational studies of the Coma 
cluster (Gregg & West 1998 and Bernstein et al. 1995). 
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Central Matter Distributions in Rich Clusters of Galaxies from z~0 to z~0.5 



If this model is confirmed, this would imply that clus- 
ters of galaxies mass profiles are well traced by the bright 
galaxy luminosity density profiles. 
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